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Abstract. We introduce and analyze a spectral vanishing viscosity approx- 
imation of periodic fractional conservation laws. The fractional part of these 
equations can be a fractional Laplacian or other non-local operators that are 
generators of pure jump Levy processes. To accommodate for shock solutions, 
we first extend to the periodic setting the Kruzkov-Alibaud entropy formu- 
lation and prove well-posedness. Then we introduce the numerical method, 
which is a non-linear Fourier Galerkin method with an additional spectral vis- 
cosity term. This type of approximation was first introduced by Tadmor for 
pure conservation laws. We prove that this non-monotone method converges 
to the entropy solution of the problem, that it retains the spectral accuracy 
of the Fourier method, and that it diagonalizes the fractional term reducing 
dramatically the computational cost induced by this term. We also derive a 
robust Z^-error estimate, and provide numerical experiments for the fractional 
Burgers' equation. 



1. Introduction 

In this paper we are concerned with a spectral vanishing viscosity (henceforth 
SVV) approximation for periodic solutions of non-local or fractional conservation 
laws of the form 

d t u + d x - f(u) = -(-A) x / 2 u, (x,t)eD T 
u(x,0) — uq(x), xgA, 



for A G (0, 2), or more generally, for 
(1.1) 



d t u + d x -f(u) = £"[«], (x,t)€D T 
u(x,0) = uq(x), x G A, 

where D T = A x (0,T) and A = (0,27r) rf , and £»[■] is a non-local (Levy type) 
operator defined as 



(1.2) m<P{-)]{x) = / <f>(x + z) - <j>{x) - z ■ d x <f>(x) 1 N<1 dfi{z), 

J\z\>0 

where l(-) is the indicator function. Throughout the paper we assume that 
(A.l) / = (/i, . . . , f d ) with fj G C S (R) for all j = 1, . . . , d (s to be defined); 

(A.2) , > is a Radon measure such that / \z\ 2 A 1 d M (z) < oo; 

J\z\>0 

(A.3) u Q G L°°{A) n BV(A), u is A-periodic. 
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Here and in the rest of the paper, a A 6 = min(a, 6), 
d d 
dt = dt ' dj = dx~ 9x = " ' ' ' dd '' 

Integro-PDEs like (|l.ip typically model anomalous convection-diffusion phenom- 
ena. When /i = n\ is dchncd by 

(1.3) dTT X {z) = j-^dz, c x > and A G (0, 2), 

then £ = — (— A) A / 2 and the equation finds applications in e.g. over-driven det- 
onation in gases [14] and anomalous diffusion in semiconductor growth [38 . Ap- 
plications in dislocation dynamics, hydrodynamics and molecular biology can also 
be found, see e.g. the references in [HUH]- Many more applications can be found 
if asymmetric measures /i are allowed. For example, we cover the (linear) option 
pricing equations for all Levy models used in mathematical finance [151 134] , if 

(1.4) dfi(z) = g(z)dir\(z), 

for some possibly asymmetric locally Lipschitz continuous function <?(•). An exam- 
ple is the one-dimensional (d — 1) CGMY model where 

' Ce- G ^ for z > 0, 

Ce- Mlzl for z < 0, 



g{z) 



for positive constants C, G, M (and Y = A). In general the non-local operator C is 
the generator of a pure jump Levy process, and conversely, any Levy process will 
have generator like C when (|A.2[) is satisfied. We refer to the books [3] [131 132] 
for more information about Levy processes and their many applications. The most 
general Levy measures for which the results of this paper applies, are Levy measures 
\i that can be decomposed as 

(1.5) n = Li s +n n , 
where 

(1.6) {J-s,Hn > 0, /i s is symmetric and / \z\ A 1 d/i n (z) < oo. 

J\z\>0 

See Section [5] for statements of results and remarks. This class possibly includes all 
Levy measures, but we have so far not found a proof of this. At least it includes 
all the Levy measures found in finance, see Remark 18.31 and also many singular 
measures like e.g. delta-measures. 

It is important to note that non-linear equations like (|1.1|) do not admit classical 
solutions in general, and that shock discontinuities can develop even from regular 
initial conditions. This is well known for pure conservations laws (where C — 0), 
see e.g. For fractional conservation laws where £ = — (— A) A / 2 , it is shown in 
recent works that solutions are smooth for A G [1,2) [HI HH HE] • However, when 
A G (0, 1), the fractional diffusion is too weak to prevent shock discontinuities from 
forming, see [1] [TO] [26]. In some cases however, these shocks are smoothed out 
over time [5]. When shocks form, weak solutions become non- unique and entropy 
conditions are needed to select the physically correct solution - the entropy solu- 
tion. The well known Kruzkov entropy solution theory for conservation laws was 
extended to fractional conservation laws in [T]. This extension relies on new ideas 
for the fractional term and is strongly influenced by the viscosity solution theory for 
fractional Hamilton-Jacobi-Bellman equations. Extensions of the Kruzkov- Alibaud 
theory to general Levy operators and even non-linear fractional terms can be found 
in [HES]. 
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In this paper we deal with a SVV (spectral vanishing viscosity) approximation 
of A-periodic entropy solutions of (|1.1[) . The method is a Fourier Galerkin method 
with an additional spectral viscosity term. Because of the formation of shocks 
in the solutions of (jl.ll) . it is very difficult to devise a convergent and spectrally 
accurate numerical approximation of this equation. This has to do with the fact 
that Fourier spectral methods support spurious Gibbs oscillations, and thus fails 
to converge strongly toward discontinuous solutions. It is well known that such 
methods need to be augmented by some kind of vanishing viscosity in order to 
achieve convergence. But the standard vanishing viscosity method is not spectrally 
accurate. To overcome these problems, we use the SVV approximation developed 
by Tadmor in [35], cf. also Q31 1301 1331 GIH] and the books g] [6]. To suppress 
spurious oscillations without sacrificing the overall spectral accuracy of the method, 
Tadmor adds a modified viscosity term, which in Fourier space only affects high 
frequencies. There are two parameters involved in this approach, the coefficient of 
the viscosity term e and the size m of the viscosity free spectrum. Spectral accuracy 
and convergence toward the unique, possibly discontinuous, entropy solution, then 
follows by imposing appropriate conditions on e and m. We also like to mention 
another important feature of the method. In all cases, it diagonalizes the fractional 
term and hence reduces dramatically the computational cost induced by this term. 
In our rather naive implementation for the fractional Burgers' equation, the SVV 
method turned out to be orders of magnitude faster than a Discontinuous Galerkin 
approximation of the same equation where the fractional term gives full matrices. 

When equation (ll.lj) is linear, f(u) = u, or when it is local C = 0, there is a vast 
literature on numerical methods and analysis, some methods and many references 
can be found e.g. in [4] [TSJ [23] [34]. In the general case however, there is not much 
work on numerical methods, we only know of the papers [131 [TTJ [2T]. Difference 
methods are introduced in [17) for equation and in |21] for an equation sim- 

ilar to (|l.ip from radiation hydrodynamics. In [17) . the first general convergence 
result for monotone schemes is obtained. Finally, in |13j . a Discontinuous Galerkin 
approximation of (|1.1[) is analyzed and a Kuznetsov type of theory is established 
and used to derive error estimates. A periodic extension of this theory will be used 
to find error estimates in this paper. 

Throughout the paper we will use the following additional notation. A sub- 
script p indicates A-periodicity in the space variables (i.e. in or Cf°). Here 
A-periodic means 27r-periodic in each coordinate direction. As a generic constant 
we use C. Note that the value of C may change from line to line and expression to 
expression. We also need notation for high order derivatives and their norms. Let 
a = (ax, . . . , ay) be a multi index, then 

d:=d^d^---dT, «;= U (4 and ||«, = 

\a\— s \a\—s 

Remember that aj > 0, |a| =«!+••• + ay, and that x a — a;" 1 • • • x^ d for any 
x e R d . 

The rest of this paper is organized as follows. In Section [2] we introduce an 
entropy formulation for periodic solutions (|l.ip . and give a /^-contraction and 
uniqueness result. In the same section we introduce the classical vanishing viscosity 
approximation of (11.11) . and show convergence towards (11.11) with optimal L 1 error 
estimate. As a corollary we get existence for (|1.1[) . The proofs rely on the Kruzkov's 
doubling of variables device [2 [27] and Kuznetsov type of arguments [13l [28] , and 
are given in the Appendix. The SVV approximation of (jl.ip is introduced in Section 
[3] and we show that it is spectrally accurate and that it diagonalizes the non-local 
operator. In sections HHS] we assume that the measure /x is symmetric. In Section^] 
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we prove an energy estimate for the SVV method. Along with results from [TT] , this 
allows us to control the so-called "truncation error", the spectral projection error 
coming from the non-linear term. In Scction[5]we prove a priori L°°, BV , and time 
regularity estimates for the SVV method, and obtain compactness. In Section [6] we 
prove that the SVV method converge to the classical vanishing viscosity method 
from Section [SJ Combined with the results of that section, it follows that the SVV 
method converges to the entropy solution of (jl.lj) . In the process, we also prove the 
optimal I^-rate of convergence for our SVV approximation. We solve numerically, 
using our SVV method, the the fractional Burgers' equation in Section [7] Finally, 
in Section [3] we extend the results in the previous sections to allow for asymmetric 
measures /i. 

2. Entropy formulation for periodic solutions 

In this section we introduce an entropy formulation for A-periodic solutions of 
the initial value problem (jTTTJ) . To this end, we write the operator £**[•] as 

m ( j>]=L^]+^' r [4>]-ll-d^ 

where 

C-r ["MOKz) = / <f>( x + z )- 4>(x) ~ z ■ d x 4>{x) 1m<! dn(z), 

J\z\<r 

C^(-)](x) = / <Kx + z)- <t>(x) dn(z), 

J\z\>r 

7£ = / z d[i{z). 

Jr<\z\<l 

If r > 1, we take 7^ = 0. The adjoint of takes the form 
where 

C* r >»[<j>(-)](x) = f cp{x~z)~ <f>{x) + z ■ d x 4>(x) 1| ZK1 dfi(z), 

J \z\ <r 

£** r [0(.)](») = f 4>(x-z)-<t>(x)dii(z). 

J \z\>r 

We also let r), 77', and q denote the functions 

j](u, k) = \u- k\, r/(u, k) = sgn(w - k), qj(u, k) = rf{u,k) (fj{u) - fj(k)). 
We now define the solution concept we will use in this paper. 

Definition 2.1. (Periodic entropy solutions) A function u is a periodic entropy 
solution of the initial value problem (jl.ll) provided that 

1) u e C([0, T]; Z£°(K d )); 

ii) for all k E M., all r > 0, and all nonnegative test functions if £ C£°(R d x 

(o,m 

// i](u, k) d t <p + q(u, k) ■ d x ip 
(2.1) JJdt 

+ n{u, k) C* r ^[ip] + r)'(u, k) C^ r [u] tp + n{u, k) 7^ • d x ip dxdt > 0; 
Hi) esslim t ^ ||"(-,i) - mo(')IIli(A) =0. 
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Remark 2.1. In the entropy inequality (j2.ip it is easy to see that all the terms, 
except possibly the £ A1 ' r -term, are well defined and A-periodic in view of i). The 
problem with the £ M ' r -term is that we integrate a Lebesgue measurable function 
w.r.t. a Radon measure /i. But the term is still well defined because the integrand 
of C^ ,r [u] is measurable w.r.t. the product measure dn(z)dxdt. This is true because 
the integrand is the dfj,(z)dxdt-a,.e. limit of continuous functions, a fact which read- 
ily follows from the fact that u is the da;di-a.e. limit of smooth functions. By i), 
(A.3), and Fubini, we then find that C' r [u] G C([0, T]; L™(R d )). 

We now state the following central result: 

Theorem 2.2. (L 1 -contraction) Let u and v be two entropy solutions of the initial 
value problem (jl.ip with initial data uq and Vq. Then, for a.e. t G (0,T) ; 

(2.2) \\u(;t)-v{;t)\\ L t (A) < \\u -v \\lHa)- 

The proof will be given in Appendix|21 Uniqueness for periodic entropy solutions 
of (jl.ll) immediately follows by setting uq — vq. 

Corollary 2.3. (Uniqueness) There is at most one entropy solution of 
We now consider the vanishing viscosity approximation of 

^ ' ' \ u e (x, 0) = Uq(x) x g A. 

In this paper we always assume that this problem admits a unique classical solution 
u c . This is of course true, but a proof lays outside the scope of this paper. Remark 
!2.6l below provides some ideas on how to prove this result. We now give an estimate 
on the rate of convergence of u £ toward the entropy solution u of (jl.ip . 

Theorem 2.4 (Convergence rate I). Let u be the periodic entropy solution of (jl.ip . 
and u e be a smooth solution of (|2.3p . Then, 

(2.4) |K,*)-u £ (-,t)|| iI(A) < Cy/l. 

The proof is given in Appendix [Bj This result generalizes to periodic fractional 
conservation laws Kuznetsov's well known result for scalar conservation laws |28j . 
As a by-product of the well-posedness of (|2.3p and Theorcm l2.41 we have the exis- 
tence of entropy solutions of (jl.ip . 

Corollary 2.5. (Existence) There exists an entropy solution of (jl.ip . 

Remark 2.6. Uniqueness of solutions of (j2.3p can be proved using an entropy for- 
mulation (see the start of Appendix [Bj and a standard adaptation of the proof of 
Theorem 12.21 incorporating ideas of Carrillo [7] to handle the Laplace term. Ex- 
istence of an entropy solution can be proven e.g. by appropriately modifying our 
spectral approximation, compactness, and convergence analysis, see the following 
sections. The solution of (j2.3D will also be smooth. To see this, note that the 
principal term in (12 .3[) is the eA-term while the /^-term is of lower order, and 
hence regularity proofs for viscous conservation laws f (|2.3p with \i = and e > 0) 
should still work after some modifications. We refer to e.g. [21] for regularity of 
viscous conservation laws, and note that the modifications typically consist of using 
interpolation inequalities for the Z^-term, see e.g. Lemma 2.2.1 in [22 . 

3. The spectral vanishing viscosity method 

We introduce a Fourier spectral method for the d-periodic initial value problem 
(jl.ip . The approximate solutions will be A^-trigonometric polynomials, 

U N ( X ,t) = ^2 "{W^l 
\t\<N 
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which solve the semi-discrete spectral vanishing viscosity (SVV) approximation 

d 

(3.1) d t u N + d x ■ P N f(u N ) = C^[u N ] + e N E d ]kQ J N * u n 
with 

(3.2) u N {x, 0) = P N u {x), 
where the Fourier projection P/v is defined as 

p N ^x) =Y,h s ' x for h = Try [ e ^' x dx - 

The (spectral) vanishing viscosity term has the following three ingredients: 
(A. 4) a vanishing viscosity amplitude cn ~ N~ e with < 9 < 1; 
(A. 5) a viscosity-free spectrum m^v ~ N% (log N)~% ; 
(A. 6) a family of viscosity kernels 

N 



;N {x,t)= ^ QJ'*(t)2^ e ' 

p— m.v ?I=P 



satisfying 



— Q^' fc is monotonically p- increasing, 

— Q^ fc spherically symmetric, Q^ k = Qy k for all |£| = p, 

— \Qy k - 8 jk \ < Cm 2 N p- 2 for all p>m N . 

Such kernels can be conveniently implemented in Fourier space, 

d N I d \ 

E * u * = - E E 3f (*) & e* 

j,fc=l |{| =mj v \J,fe = l / 

Combined with one's favorite ODE solver (e.g. Euler, Runge-Kutta, etc.), (|3.1j) and 
(|3.2|) give a fully discrete numerical approximation method for (II. ip . 

With left-hand sides set to zero (fi = and £jv = 0), (|3.1[) becomes the stan- 
dard Fourier approximation of (|1.1[) . It is well known that this approximation 
is spectrally accurate but, as opposed to the equation, it lacks entropy dissipa- 
tion. The approximation supports spurious Gibbs oscillations which prevent strong 
convergence toward solutions containing shock discontinuities. If the vC^-term is 
present in the equations, shock solutions are still possible in some situations [5], 
and the problem of the Gibbs oscillations remains. In order to suppress such os- 
cillations without sacrificing the overall spectral accuracy of the method, we have 
followed Tadmor [35) and added a vanishing spectral viscosity term to the scheme, 

An important feature of Fourier method (|3.1[) is that it diagonalizes, and hence 
localizes, the non-local operator This leads to dramatically reduced compu- 

tational cost for this term. Indeed, 



(3. 3) C[un]= E G"(0fi*(t)e**, 

\£\<N 

where 



(3.4) G M (0 = / e**-l-i£-zl w<1 d/i(aO. 

|z|>0 
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Furthermore, when the measure /x is symmetric, 

(3.5) K~ B ) = for a11 Borel sets B eR d \ {0}, 

the weights (|3.4p are all real and non-positive. This follows since the imaginary 
part of the integrand is odd and the real part is even and non-positive {e^' z = 
cos(£ • z) + i sin(£ • z)). Finally, we stress that the approximation of the non-local 
operator (jl.2l) is spectrally accurate since, by Taylor's formula, 

\\C^U N (;t)}-C^ U (;t)}\\ L2{A) 

< C ^SUp||a,<9 fc (wAr - «)(',*) || £2 ( A ) + ||(ujv - w)(-,i)||£=(A)J- 



Now we define 

(t) = 

and note that 



^fc-Qf(t) |£|>mjv, 



ISI<JV 



(3.6) Au N (;t) = E djdkQtfi-it) * U N (;t) + E d 3 d k R% k (;t) 
3,k=l ],k=l 

To conclude this section, we recall that by Lemma 3.1 and Corollary 3.2 of [TT], the 
spectral vanishing viscosity term is an L p -bounded perturbation of the standard 
vanishing viscosity cnAun: 

Lemma 3.1. For < r < s < 2, 

(3.7) 



E m '"v 



< Cm% {\ogN) d . 



£i(A) 



Moreover, if cn < Ce N m N (\ogN) d < C, then for all p >!,<£>€ L P (A), 



j,k=l 



< CN IMIlp(A)- 



Lp {A) 



4. Spectrally small truncation error for symmetric // 



(3.8) 



In this section we assume that the measure /x is symmetric, cf. (|3.5p . In the SVV 
approximation (|3.1[> . the convection term d x ■ f(u) is replaced by d x ■ -P/v/(wjv) 
which leads to the (truncation) term error 

d x -(I-P N )f(u N ). 

We will now show that this error is spectrally small due to the presence of the 
spectral vanishing viscosity term. 

Let us start by noting that a straightforward estimate leads to 

iiw-p^)/(^)II£= ( a) = (e E iriV^)(0i 2 )"< ''^ + ^; )llL2(A) 

for all multi- indices a, (3. Note that there is no divergence in this estimate, so d"f 
is a vector. By Theorem 7.1 in |11) . there is a constant 1C S such that 



(4.1) \\d s J{u N )\W{A) < K s \\d s x u N \\ L 2 ( A) for /C s <C7El/l 



i life — i 

C k \\ u N\\ L ^^ A ) 



fc=l 
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and s = 1,2,..., where \ f\ c * = ||^/0)IU-(n N ) and £l N = {u : \u\ < ||ujv|| L oo (a) }. 
This inequality is a type of Gagliardo-Nirenberg-Moser estimate, and similar results 
can be found in page 22 in Taylor [3 7) . By these two inequalities we can conclude 
that, for all < r < s, 

(4.2) \\d r x (I-P N )f(u N )\\ LHA) <^-\\a^ UN \\ LHA) , 

Inequality (|4.2I) states that the r-derivative of the truncation error decays as 
rapidly as the s-smoothness of un permits. Of course the s-derivatives of an ar- 
bitrary iV-trigonometric polynomial un may grow as fast as N s , in which case 
nothing is gained from (|4.2p . However, if ujy is solves our VVS approximation 
(|3.1[) . we can have the better bound ej/ in L 2 . This will be a consequence of the 
following energy estimate: 

Theorem 4.1. Consider the SVV approximation (|3.ip with o,nd itin such that 

( 8d«/C s+1 
(A.7) 6N > "^V - ' 

[e N m 2 N (\ogN) d <C. 

Then there is a constant B s (proportional to n s k=l 1C s for s > 1 and to ||uat||l°° for 
s = 0) such that 



(4.3) 



eJ«u*(-,t)IU*(A)+e5r ("EE ^(Olff f Mr)\ 2 dr 

+ 6^||^ +1 UA r|| i2(z5T) <6 s + 3e^||c)> A r(.,0)|| i2(A) . 



Remember that in this section [i is symmetric and hence G M is real and non- 
positive. Now if 

(A. 8) |/|c s < 00 for sufficiently large s, cf. (|4.7I) below, and 
(A. 9) u is such that \\d^u N (-, 0)|| L 2( A ) < C, 



then Theorem 14. II implies that 



ll^>Jv(-,i)llL2(A) < Ce N s and ||^ +i UAr|| L 2 (Z5T) < Ce^ 
Taking into account (|4.2[) , we then find that 

(4.4) ||^(J-P w )/( Ujv (-,i))|| i2(A) <CB a N~", s r = s(l-9)-r, 

(4.5) ||^(/ - P w )/(wat)|| l2(Z5t) <C6 s 7V-( s "+i), Vs > 1. 

We can now turn these inequalities into spectral decay estimates in the uniform 
norm using the Sobolev inequality (cf. Theorem 6, Chapter 5, in |20j ) 

lia^lU- <c\\d r x +l * ]+1 <p\\ L 2. 

For example, inequality (|4.5[) becomes 

(4.6) \\d r x (I ~ P N )f(u N )\\ L ~ (DT) <CB s N- s >-+W +1 --i <CB S N- S ^ +1 . 

Note that the polynomial decay rate in (|4.6p can be made as large as the C s - 
smoothness of /(•) permits. Taking r — 2, we can find the following result. 
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Theorem 4.2. If f g C s with 
(4.7) 

then 



s > 



4+[#] 



(4-8) 115,(7 - P N ).f(uN)\\L^(D T ) + - Pn)I(un)\\ L oo {Dt) < 



CB S 
N ' 

The smoothness requirement (|4.7p will be sufficient for all the estimates derived 
throughout the paper. 

Proof of Theorem \4- 1\ For sake of brevity, we will write || • || instead of || • ||l 2 (A)- 
With (|3.6p in mind, we rewrite the SVV approximation (|3.ip in the two equivalent 
forms 



(4.10) 



(4.9) d t u N + d x ■ P N f(u N ) - C[u N } - cwAmjv = -e N E djd k R% k * u N , 

3,k=l 

d t u N + d x ■ f{u N ) - C^[u N ] - e N Au N 

d 

= -ejv E d J d kR J N * u N + d x ■ {I - Pjv)/(ujv). 
Since G At (^) < (/x is symmetric) and un{x) and /^[uat] are real, 
f da: = ^ G"(0l«f(*)l 2 < °> 

•* A \£\<N 

and hence spatial integration of (|4. 10[) against ujv yields 

Id 

2d* 1 



--!!'^|| 2 - E G "(0|fi6(*)| 2 +^11^11 



< £Ar||uAr|| 



l€l<JV 
d 



E +E ii^iiik/- 

j",fe=l 3=1 



Using (I3.8P with p = 2 for the first term on the right and (I4.2p with (r, s) = (0, 1) 
for the second term, we find that 



Id 

2d7 l|uw| 



|£|<iV 



iV 



with ejv < Ce N m 2 N (logN) d < C. Hence (|43)) follows for s = since by (|A.7p . 



AT 



> 



and c w ||wAr|| 2 < C ||ujv|li«.( A ) = #o- 

The general case follows by induction on s. Spatial integration of f)4.9[) against 
c^ q ujv for some multi- index a yields 



Id 

2dt 



(4.11) 



ii^h 2 - e G»{m a \ 2 \ut(t)\ 2 +tN\\d«d x u N \ 

\i\<n 

d 



< e N \\d x u N \\ 



E d d k R% k * 9« 



l^9, UJ v||||ai Q i- 1 a,-p 7V /(u w )||. 
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After having used (13.81) and Young's inequality to bound the first and second term 
on the right hand side, we hnd that 

^ t W d >Nf- E GnO\C\ 2 \H(t)\ 2 + e f\\d«d x u N \\ 2 

<C||9> JV || 2 + -L\\dl"\p N f( UN )\\ 2 . 

Now we sum over all \a\ = s to hnd that 
1 d 

-w>N\r - 

(4.12) |a|=«|£|<2V 



2^n^ w |i 2 - E E G»m a \ 2 \udt)\ 2 + ^\\d s x +1 u N \ 



<C\\d* UN \\ 2 + £-\\d°P N f(u N )\\ 2 . 

By g3J and gj), 

\\d s x P N f(u N )\\ < \\d s J(u N )\\ + \\d s x (I - P N )f(u N )\\ 

<MaM + ^||^ +1 UJV ||, 

and hence by inequality (|4.12p we see that 
(4-13) 

H\\9>n\\ 2 - e E G^om 2 \hm 2 +( e f-^^)\\d s x +1 u N \\ 2 

\a\=s \(\<N ^ N ' 

< (c+££) \\d°u N \\ 2 <^||9> W || 2 , 
\ ejv / cat 

where the last inequality holds for N big enough. By (|A.7[) and integration in time, 
we then hnd that 

(4.14) 

hd*u N (;t)f - E E G "(0\C\ 2 f\udr)\ 2 Ar+"-fm +1 u N \\l 2{DT) 

\a\=8\£\<N 

9d s lC 2 1 

< s -\\d s x u N \\ 2 2{DT) + -\\d s x u N (;0)f- 

At this point (|4.3p follows by the induction assumption on s since 

\\d s x u N \\% (DT) < CBUe^ 211 ^. 
The proof is now complete. □ 



5. A PRIORI ESTIMATES AND COMPACTNESS 

In this section we prove uniform 

L°°(D T ), L°°(0,T;BV(A)), and C° ^([0, T];L 1 (A)) 

bounds on the solutions {un ■ N £ N} of the SVV approximation (|3.ip . As a 
consequence we obtain compactness in L . 
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5.1. Regularity in space. 

Lemma 5.1. (£°°-stability) Let (A.1)-(A.9) and (|3.5p hold anduN be the solution 
of the SVV approximation (|3.1I) . Then for t < ChxN, 

\\u N {-,t)\\L°°(A) < C||ujv(-,0)||i«( A ). 

Proof. For sake of brevity, we write just || • ||oo instead of || • ||l°°(a)- First we note 
that, for any smooth convex function rj{-) with derivative /?'(■), we have that 



(5.1) 



V '(u N )£^[u N }<^[r,(u N )}. 



This is a consequence of the inequality jf(b)(a — b) < 77(a) — rj(b) which holds for 
all smooth convex functions rj(-). Moreover, 

(5.2) f C[ri(u N (;t))](x) dx = 0. 

J A 

To see this note that 

\r](u N {x + z)) - rj(u N (x)) + z ■ d x rj{u N (x)) l|d<i I dfi(z) dx 



A J \z\>0 



< \\d 2 x v(u N )\\ 



\z\ 2 d/i{z) + ||??(ujv)|| 



z\<l 



d[i(z) < oo, 



\z\>l 



since un is smooth and periodic. By Fubini we then find that 
£"fo(utf(-,t))](sO dx 

/ t]{u N (x + z)) - r](u N (x)) + z ■ d x rj(u N (x)) l\ z \<i dxd[x(z). 

\z\>0 J A 

By A-periodicity of , (|5.2I) now follows since 



rj(u N (x + z))dx = I Tj(u N (x)) 
J A 



dx 



for every z, and 

d Xi r]{u N (x' , Xi)) dx'dxi 

t](un(x', 2tt)) dx' — 



(0,2tt) 



r](u N (x',0))dx' = 0. 



Let us now integrate (I4.10P against the function p u p N 1 (with p even) , and use 



'(0,27r) d - 1 

r integra' 

(|5.ip and (I5.2[) to get rid of the non-local operator £ M [-]. We then find that 

Pll«iv(-,*)Hi7(A) fo\\ U N(-,t)\\LP(A) = ^II m v(-,0IIlp(A) = J^U P N ' 1 (x,t)dtU N (x,t)dc 

<p u p N 1 (x,t) e N Y] djdkR% k *u N (x,t) + d x - (I - P N )f(u N (x,i)) dx 
which by the Holder inequality (with p and q = r^j) is less than or equal to 



EAT 



"Lp-i (A) 

djd k R% k * u N (-,t) 

j,k=l 



LP (A) 



+ \\d x -(I-P N )f(uN(;t))\\ LP{A) 
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^_e_ — H^ll^p 1 , we may divide both sides by p\\uN(-,t)\\ p LP ^ and 



Since H^" 1 

send p — > oo to discover that 
d 



dt 



||«jv(-,*)||oo < e 



iV 



^ djd k R% k * u N (-,t) 

3,k=l 



+ \\d x -(I-P N )f(u N (;t))\\ c 



By ([3~g]l . ([575)1. the definitions of B s , /C s and cjv, and (A. 7), it follows that 

R • (J ~ PA)/(^(-,i))||oo < § < ^ J! /C s < Slujvlli, 

fc=i 



f/V 



djd k R% k (-,t) * ujv(-,t) 



< C\\u N \\ 



and hence 



d C 

-7rll"Ar(-,i)||oo < cjv||wjv(-)*)||oo + tt ||wjv(*, 
dt Jy 



Letting y(t) 
we find that 



i)||ooj and multiplying by the integrating factor e Cjvt , 



which implies that 

y(t) < y(0) ll 

Going back to ||itjv(-, £)||ooj we can conclude that 

||wAr(-,t)||oo < e C «*||-U A r(-,0)| 



4r-l 



□ 



where the last factor is bounded for £ < ClniV for some C. 

We also have the following result: 

Lemma 5.2. (BV-stability) Let (A.1)-(A.9) and (|3.5p ZioW, and un be the solu- 
tion of the SVV approximation (|3.ip . TTien 

||«w(-,r)|| flV(A) < e c " T (\\u N (; 0)||w(A) + C AT S ^ 

im£/i cat = e N m 2 N (logiV) d < C and S2 = s(l - 6) - 2 > 0. 
Proof. Spatial differentiation of (|4.10[) yields 

d t diU N + d x ■ (/' \u N )diU N ) - C^ldiUN] - e N Ad l u N 

d 

= did x ■ (I - Pn)J(un) + eat ^ 'hi), J, ,J y * d l u N . 

If we integrate this expression against sgn e (9jUjv), where sgn g (-) is a smooth ap- 
proximation of the sign function, we can get rid of the non-local operator as 
in the proof of Lemma 15.11 If we also use (|3.8[) with p — 1 and take the limit as 
g — > 0, a standard computations reveal that 



— \\diU N (;t)\\ L i^ < C\\did x ■ {I - P N )f(u N )\\ L x^) +cjv||9iWjv(-,*)|Ui( A )- 



PERIODIC FRACTIONAL CONSERVATION LAWS 



13 



Since ||ujv( - ) £)||.BVTA) < J2i=i H^iUzvC ^)IIl 1 (A)j we integrate this inequality in time 
to see that 

\\uN(-,t)\\ BV(A) < e CNt I \\u N (; 0)|| W(A) + C ||^(/ - P N )f(u N )\\mD T ) ) ■ 



But by (1441) . 

||^(/ - P N )f(uN)\\mD T ) < C \\dl(I - P N )f(u N )\\ LHDT) <CB S Vtn~ s \ 
and the proof is complete. □ 
5.2. Regularity in time. 



Lemma 5.3. (Regularity in time) Let (A.1)-(A.9) and (|3.5p hold and, un be the 
solution of the SVV approximation (|3.ip . Then 



\\u N {-,ti) - u N (-,t 2 )\\ L i( A ) < C - h\. 

Proof. Let u N {-,t) — Wjv(-,i) * w e (-) for an approximate unit co e (cf. the proof of 
Theorem 12. 2[) . By the triangle inequality we see that 



(5.3) 



\\u N (-,h) - U N (-,t 2 )\\ L i(A) < ||Wjv(",*l) ~ Mtf (' , *1 ) || L 1 (A) 



+ \\ u N(-,h) - U e N (-,h)\\n(A) + \\u N {-,t 2 ) ™ MlV (' , h) \\ V (A) ■ 

The first and the third term on the right-hand side of (|5.3[) are bounded by e|u|sy: 



\\u N {-,t) - u^(-,i)|| L i (A) 



uj e (y - x)(u N (x,t) - u N (y,t) ) dy 



d.r 



< 



uj e (s) un{x, t) — un{s + X, t) 



ds dx 



I A JR d 

<^fd\u{-,t)\ BV(K) I \s\u) t {s)ds 

< Vde \u(-,t)\ BV{A y 
Let us estimate the second term. By Taylor's formula with integral remainder, 
\\u N (-,h) -u e N (-,t 2 )\\ L i(A) 



< \ti-h\ 



Idtw^Mi + r{t 2 - h))\ drdx. 



A Jo 



We now derive a bound for ||<9tiijv||.E, 1 (and hence also for ||9tu^||ii) by using the 
SVV approximation (|3.1[) itself. To this end, we take the convolution product of 
both sides of (14.91) with uj e to obtain 



\\dtu%\\ L i {A) < \\d x ■ P N f{u N ) * w c ||li(A) + W^Iun] * w e || L i (A) 



+ e N \\Au N * w e |Ui( A ) + e N 



^2 d 3 d k R% k * un J * w e 



Li(A) 



= h+h + h + h- 
By the triangle inequality and Young's inequality for convolutions, 
h = \\0 X ■ PNf(uN) * W £ || £ i (A) 

< || ^a; • f(u N ) * W e || L i( A ) + \\d x ■ (I - P N )f(u N ) * W £ ||li(A) 

< ||9x • /(wat)IUi(a) + \\dx ■ (7- Pjv)/(wjv)||i,i(a)- 
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Therefore, by the regularity of / and un ((A. 8), Lemmas 15.1115. 2[) and (|4.8|) . we 
find that 

Il<c(\u(;t)\ BV(A) + ^y 

For the term containing the non-local operator we write 



h < 



un{x + z) — un(x) — z ■ 8 x un(x) dfj,(z) uj e (x — y) dy 



z <1 



da; 



Un{x + z) — un(x) dfi(z) u) e (x — y) dy 



z >1 



dx 



The second term on the right hand side of the inequality above is easily seen to 
be bounded by C\\un{-, fyWt 1 > while Taylor's formula with integral reminder and 
integration by parts reveals that the first term is bounded by 
-i 

(1 - t) \z\ 2 \d x u N (x,t)\ \d x u e (x - y)\ drd/j(z)dydx 



z\<l Jo 



< Ce- 1 \u 



BV(A)- 



For the Laplace term we have 

^3 < \\d x u * a x w e || L i( A ) < e^ 1 \u\ B v(A), 
and finally, using Young's inequality for convolutions and 

d 

^2 djd k R J N k * u N j ■> ^, 



h — ejv 
To sum up we have 



< C\\u N \\ L i {A) . 



Li(A) 



\\dtu e N \\ L i {A) < \\d t u N \\ L i {A) < C ^1 + -^j , 
and inequality (|5.3p and the above estimates then implies that 
\\u N (-M)-u N (-,t 2 )\\L^A) <CU+|ti-i2| + 
Take e — yj\t\ — t 2 \ and the proof is complete. 



□ 



5.3. Compactness. Thanks to the space/time a priori estimates in Lemmas 15.11 
- 15.31 and a Helly like compactness theorem, cf. Theorem A. 8 in [23], the family 
{it at : N £ N} of solutions of the SVV approximation (|3.ip is compact. 

Theorem 5.4 (Compactness). Let (A.1)-(A.9) and (|3.5p hold, andu^ be the solu- 
tion of the SVV approximation (|3.1[) . Then there exists a subsequence un converg- 
ing in C([0,T]; L X (A)) to a limit u £ C([0, T]; L 1 (A))nL oo (_DT)nL oo (0, T; BV(A)). 

6. Convergence and error estimate 

The solution v tN of the vanishing viscosity method (I2.3[) converges to the unique 
entropy solution u of (jl.ll) . and by Theorem l2.41 

II «(•>*) - v tN {-,t)\\ L i {A) < CVfjv. 

In this section we prove a similar error estimate between v eN and the SVV approx- 
imation Mat. 
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Theorem 6.1. Let (A.1)-(A.9) and (|3.5[) hold, un be the solution of the SVV 
method (|2.3D . and v tN be the solution of (|3.1j) . Then 

\\u N (-,T) - w £jv (-,T)|| L i( A) < C y/ejf. 

A direct consequence of Theorems 12.41 and 16.11 is the following convergence and 
error estimate for the SVV method. 

Corollary 6.2. (Convergence with rate) Let (A.1)-(A.9) and (|3.5[) hold, u A be 
the solution of the SVV method (|3.ip . andu be an entropy solution of (jl.ll) . Then 

\\u(-,T) - u N (-,T)\\ L i( A} < C y/etf- 

Proof of Theorem \6.1\ Since v tN is smooth, we can subtract equation (|2.3[) from 
equation (|3.ip to obtain 

d t (u N -V CN )+d x - (/(Un) - /(«£«)) - C^UN - Ve N ] ~ eN&{ u N ~ 

d 

= -*N ^2 d i R3 N * d kV>N + 9 X (I - Pn)/(un)- 

As explained in the proof of Lemma l5.21 we can integrate such an inequality against 
(a smooth approximation of) sgn(u7v — v €N ), to find that (after going to the limit) 

: ||ujv - v 6w ||li(a) 



dt 

< 



e N \\ d 3 R N \-,t)*d k u N (-,t) +\\d x -(I-P N )f(u N (-,t))\\ L i {A) . 
i,fe=l L (A) 



Li(A) 



By ([XT]) with r = s = 1, (A.4), (A.5), and Lemma 

d d 
|| £ ^/(-,t)*d fcUA r(.,t) <|| £ djBtf^t) d k U N (:t) 

< C m N (logN) d \\u N (-,t)\\ BV(A) < Ce"*, 
so we can integrate in time to obtain 

\\u N (-,t)-v tN (;t)\\ LHA) < \\d x ■ (I-P n )/(u n (;T))\\ l i {Dt) 

< c(yfa + \\d x ■ (I-P n )/(u n (;T))\\ l2{Dt \ 

By g3]), 

\\d x ■ (L ^ P N )f(u N (-,T))\\ L 2 {DT) < C KC s N-^+^ < C IC S N-^ = C 
since Si = s(l — 9) — 1 > 0, cf. (|4.7I) . The proof is now complete. □ 

7. AN APPLICATION: THE FRACTIONAL BURGERS' EQUATION 

In this section we apply the results of the previous sections to numerically solve 
the fractional (or fractal) Burgers' equation in M. d , 

(71) f d t u + uj: d j=1 d Xj u = -(-A)^u, (x,t)eD T , 

\ u(x, 0) = uq(x), x G A, 

where the fractional Laplacian term — (— A) a / 2 wjv = £* x {un} and n\ has been 
defined in (| 1 . 3|) . In this setting expression (I3.4[) becomes 

5 «-*-l-^.«l w<1 



\z\>a 



\d+\ ■ 
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(a) A = 1.6 



(b) A = 1.1 



(c) A = 0.6 



(d) A = 0.1 



Figure 1. Solutions of system (|7~3T) with N = 256 and T = 0.5. 
The piecewise constant initial datum is uq(x) — sgn(7r — x). 

with c x = AT(^) (2tt5+ a T(l - f )) , cf. 19 j. We have the following result: 
Proposition 7.1. 

- C x |£| A for d = 1, 

(7.2) G^(£)={ „ _ A 



dS v for d > I, 



2/1 = 1 



where C* A = 2c\ A -1 x~ x sin x dx > and f |=1 dS a = 2 7 r d / 2 r- 1 ( 



-If d-\ 

2' 



The proof is given at the end of this section. In the above result and in the 
following, dS y will denote the surface area measure of the unit sphere \y\ = 1. 
Expression (|7.2j) is the "Fourier symbol" of the fractional Laplace operator in our 
periodic setting. When A G (0, 1), the integral 9a = J x~ x sin x dx is a generalized 
Fresnel integral |29j with value 

'Tr(l-A)' 



T(l - A) sin 



When A = 1, 9a is a Dirichlet integral [53] and has value f. For A € (1,2), the 
integral 9a has to be evaluated numerically since explicit formulas are not available. 

Remark 7.2. By Proposition 17. II there is a positive constant such that 



C n >~[u N (;t)}u N (;t) da; 



ISI<JV 



where right-hand side is a fractional Sobolev semi-norm [4] 

x; i£hM*)i 2 = iM-,*)&v W 

l«l<JV 

Simple energy estimates can then be used to show that the solutions of (|7.1[) be- 
long to if A / 2 (A), which is more regularity than what can be expected for general 
solutions of the pure Burgers' equation (fj, = 0). 
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(c) A = 0.6 (d) A = 0.1 



FIGURE 2. Solutions of system ([7l3]l with N = 256, T = 0.5, and 
ejv = 0. The piecewise constant initial datum is uq(x) — sgn(7r— x). 



We now use the SVV method (I3.ip to work out some approximate solutions of 
the fractional Burgers' equation (|7.1[) with d = 1. Hence f(u) — u 2 /2 and n = n\ 
in (|3.ip . We multiply both sides of (|3.1[) by e~ l * x , and integrate over (0, 27r) to 
obtain the following system of ODEs 

(7.3) 

— fi 4 (t) + ^ X! "p(*)^(*) + c xl^(*) + e ivi mw <iei<jv|ei 2 4(*)%(*) = o ! 

\p],\l\<N 

p+q=i 

where the Fourier coefficients satisfy the assumptions listed in Section [3] and are 
chosen as in |30j (they vary continuously between zero and one) . In our simulations 
we have used a fourth order Runge-Kutta solver for (|7.3[) . 

The results of our numerical simulations can be found in Figure [T] and [2J The 
results in Figure [T] confirm the convergence of our SVV approximation (|7.3[) for all 
all values of A e (0, 2). In Figure[5]we have solved the (|7.3p with e« = (no spectral 
vanishing viscosity). For A > 1, convergence continues to hold, while for A < 1, 
convergence fails and spurious Gibbs oscillations appear. This is consistent with the 
theoretical results for fractional conservation laws [31 [IS] : These equations admit 
smooth solutions for A > 1 (the strong diffusion case), while shock discontinuities 
may appear for A < 1 (the weak diffusion case). 



Proof of Proposition \7.1\ Let us prove the case d = 1 first. By Euler's formula, 
e i£z _ C os(^z) + i sin(£z), we find that 

f jto-l-itz f cos(gz) - 1 f sin(frQ-£z 

y w <i \z\ i+x K\<i \a i+x i\z\<i m i+a 
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Taylor expansions show that these integrals arc finite. In fact, the sin-integral is 
zero since the its integrand is odd. Integration by parts then leads to 



cos(£z) — 1 



z <i 



1+A 



dz = 2 



1 cos(^) - 1 







rl + A 



2 

-^r(cos(£z) - 1) 



1 2£ f 1 sin(£z) 



o A J 



dz 



2£ f 1 sin(£z) 



= --(cos(0 - 1) - f 

Now we consider the integral over \z\ > r. Again the imaginary part (the sine part) 
is zero, and a computation like the one we performed above reveals that 

e^ z - 1 , 2, _ _ 2£ f 00 sin(£z) 



dz= _(cos(0-l)- v 



dz. 



Note the + sign of the cosine-term! We add these two equations and find that 

2£c A [°°M&) 



A 



dz. 



The integral Jq 00 z A sin(£z) dz is finite and positive for all A g (0, 2) (cf. [TB] for 
details). Whenever £ > 0, we can use the change of variable £z — > x to deduce that 



Jo 



and thus 



2ca . A 



sini 



da; 



SIM 



da;. 



When £ < 0, we use the relation sin(— £x) = — sin(£a;) to obtain 

sin x 



G^(0 = -^IC| A 



,.A 



da;, 



and the conclusion for d — 1 follows. 

When d > 1 we use polar coordinates x = ry for r > and \y\ — 1, and we find 
that 



\z\<l 



pii-z 



kl>i 



Id+A 



dz = 



dz 



2/1=1 JO 



1 <*»(£• vr)-!^^ 



V\ = 1J1 



cos(£ • yr) — 1 

t-I+A 



dr dS,, 



Proceeding as in the d = 1 case for the r-integral with y fixed, we find that 



|y|=i 



SIM 



A 



dS 







da; 



SIM 



da;. 



'|y|=i 

By symmetry, the value of the y-integral is the same for any £. Therefore, 

A 



|y|=i 



dS,, 



/ |y • y\ x ds y = / ds y . 

J \y\=i J\y\=i 



The proof for the case d > 1 is now complete. 



□ 
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8. Extension to asymmetric measures \x 

In this section we show how to modify the arguments of the previous sections 
to obtain results for a large class of non-symmetric measures /x including all the 
Levy measures used in finance. A careful look at the previous arguments shows 
that symmetry of /x is used for the sole purpose of having a sign of the fractional 
term in the energy inequality (see (|4.14|Q in order to prove Theorems 14.11 and 14.21 
This fractional term is 

(8.i) [[ m UN ]di a u N dx-dt = J2 (o in 2 f T K(t)i 2 dt, 

JJd t ^ n Jo 

and it is non-positive when /x is symmetric. In the general case the sign of the 
fractional term (|8.1[) is unknown, but everything still works if we assume that 

fi = fi s + fjtn, 

for /x s ,/i„ satisfying (|1.6p (i.e. we assume (| 1 . 5[) and (|1.6[) ). Note that in this case, 
we may split the weights in (|3.4I) into their symmetric and non-symmetric parts, 

where G^ a (£) is again real and non-positive, and by (|1.6|) , 



f e**-l-i£-zl w< i dfi n (z) <C n (l + \S\) 

J\z\>0 V 7 



(8.2) \G^(0\ = 

■ '\z\>0 

The main result of this section is the following: 

Theorem 8.1. (Convergence with rate) Let (A.1)-(A.9), (jl.5[) and (|1.6[) hold, um 
be the solution of the SVV method (|3.ip . and u be an entropy solution of (jl.lj) . 
Then, 

\\u(-,T) - u N (-,T)\\ L i (A) < C^/e^. 

To prove this result, we have to modify the arguments of the previous sections. 
In view of the above discussion the key result to obtain is a version of Theorem 14. II 
for measures /i satisfying (ll.5[) and (ll.6[) : 

Theorem 8.2. Assume (A.1)-(A.7), (ll.5[) . (|1.6p hold, and let un be the solution 
of the SVV approximation (|3.ip . Then there exists a constant B s (proportional to 
1 + II|. =1 /C S for s > 1 and to ||ujv||l°° for s — 0, see Theorem \4-l\ l such that 

e^\\d s x u N (;t)\\ L , {A) + e s ^\\d^WhHD T ) < B s + ^ N W x u N (-,0)\\ LHA) . 

We prove this result at the end of this section. Now if we also assume that (A. 8) 
and (A. 9) hold, then it easily follows that Theorem 14.21 still holds if we replace B s 
by B s . At this point the reader may easily check that all the other results also hold 
if we everywhere replace B s by B s - and hence Theorem 18.11 follows . 

Remark 8.3. A Levy measure [i defined by 

d/i = g(z) d7r A (z), 

(see (|1.4[) ) can be written as /x = fi s + ji n where 

d^is = g(z) A g(~z)dir x and dfx n = [g(z) — g(z) A g(—z)] dTt\. 

Note that /i s ,/i n > 0, fi s is symmetric, and that \i n satisfies the integrability con- 
dition in (|1.6p if g is locally Lipschitz: Let g n (z) = g(z) — g(z) A g(— z) and note 
that g n (0) = 0, hence g n (z) = \g n (z) - 9n{0)\ < C \z\ for \z\ < 1. 

We now show how to modify the proof of Theorem 14.11 to prove Theorem 
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Proof of Theorem \8.2[ Once again we use the shorthand || • || instead of || • ||l 2 (A)j 
and rewrite the SVV approximation (|3TT]) as in and PTTU|) . Note that (ISTTj) 

and (|5.2j) holds for general measures fj,, so we find that 



/ Ciu^UN dx < / C[u%] da; = 0. 

J A J A 



Hence, spatial integration of (|4.10p against un yields 



1 d ll 



ejv \\d x u N \ 
d 



<e N \\u N \\ E djd k R% k * u N + E 1 1 dj u N \ \ || ( I - P N ) fj (u jy ) 1 1 , 

j,k=l j=l 

and the conclusion in the case s = follows exactly as in the first part of the proof 
of Theorem 14.11 

Now let s > 0, and note that by (|8.2I) and Young's inequality, 



f di a u N c^[u N \ dx= ho 2q g^(oi%wi 2 

|£|<JV 



|£|<A 



iei<jv 



If we take this into account and perform spatial integration of f|4.9[) against d 2a UM 
for some multi- index a, we find the following modified version of (I4.11[) , 



Ld, i*swir- E G^me\ 2 \udt)\ 2 + ^ L \\dZd^ N f 

\£\<N 



2dr 



< e N \\d"u N \\ 



UN 

3,k=l 



l^u^HIISH- 1 ^ • P N f(u N )\\ + ^\\d«u N \ 



As in the proof of Theorem 14. 11 we now use p. 81) and Young's inequality to bound 
the first and second term on the right hand side. The result is that 

1 d i9>ivii 2 - e c"(om a ific(*)i s 

\£\<N 

£jv " ~ ca 

Now we sum over all \a\ = s to find that 



2dr 



2 1 £ f\Kd xUN \\ 2 



< C \\d%u N \\ 2 + ±\\ d ^P N f(u N )\\ 2 + ^%\\d%u N f 



1 d 



^ t P>N\\ 2 - E E G^(ori 2 Kwi 2 + ^ii^ +i 



|a|=s \£\<N 



ejv, 
2 



2 , " I|md__j;/..._m|2 , ZL/ nnn« 



<cii9> JV r+— \\d°p N f( UN )\\ 



5>a|| 2 . 



Thanks to gT]) and 

< ic s \\d s x u N \\ 



ii os4-1 ii 
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and hence 

a| =s |£|<jV ^ N ' 

/ + 2Cl + 2dVq\ (| a ^ 2C* + 3d a Eg || 5>jv || 2) 



where the last inequality holds for iV big enough. 
To conclude, we use (|A.7|) to obtain 

\a\=8\(\<N 

, 2 Cl + 3d s /C2 ll9 .1, 



u n\\ 2 L 2(d t ) 



The proof is now complete since by induction on s, 

2 / r< R2 -(2s-1) 



ii9>^iii 2(DT ) + ? ii9>jv(-,o)r. 



□ 



Appendix A. Proof of Theorem 12.21 

Let us take ip = ip(x,y,t, s), it = and w = v(y,s). We set fc = v(y,s) in 

the entropy inequality for u(x,i), and integrate over all {y, s) £ Qt to obtain 




i]{u{x, t), v(y, s)) d t ip(x, y, t, s) 



+ q{u(x, t),v(y, s j) ■ d x i/j(x, y, t, s) 

+ r 1 (u(x,t),v(y,s))£;^[i,(;y i t,s)](x) 

+ V'(u(x, t),v{y, s)) C^ r [u{-, t)]{x) ip(x, y, t, s) 

+ i](u(x,t),v(y,s))Y p , ■ d x i/j(x,y,t,s) dxdtdyds > 0. 

In the entropy inequality for v(y, s), we set k = u{x,t) and integrate with respect 
to (a;, t) to find that 

// // l{u{x,t) 1 v{y 1 s))d s ^{x,y 1 t,s) 

JJ D T JJ D T 

+ q(u(x, t), v(y, s)) ■ d y ip(x, y, t, s) 

+ rj(u(x, t), v(y, s)) C*'» [ijj(x, ; t, s)] (y) 

- t]'(u(x, t), v(y, s)) C^' r [v{-, s)](y) ij}(x, y, t, s) 

+ r)(u(x,t),v(y, s)) 7^ • d y ip(x,y,t, s) dydsdxdt > 0. 

In the following we need the K 2d -operators 

C"' r [(/>{-, -)}(x,y) = I <j>(x + z,y + z)-4>{x,y) dfi(z), 

J\z\>r 

t*' IJ " r [(j){-, -)]{x,y) = I cp(x-z,y-z)-4>(x,y) dfj,(z). 

J\z\>r 



22 



S. CIFANI AND E.R. JAKOBSEN 



With these definitions in mind, we add the two inequalities above and change the 
order of integration to find that 



r}(u(x, t), v(y, a)) (d t + d a )i/j(x, y, t, s) 

I D T J J D T 

+ q(u(x, t), v(y, a)) ■ (d x + d y )ip(x, y, t, s) 

+ r](u(x, t), v(y, a)) C** [V>0, y, t, a)] (x) 

+ r){u{x, t), v(y, a)) [V>(x, •, t, a)] (y) 

+ rf{u{x, t), v(y, a)) &' r [u(-, t) - v(-,s)](x, y) i>(x, y, t, s) 

+ V(u(x, t), v(y, a)) 7^ • (d x + d y )ip(x, y, t, s) dw > 0. 

Here and in the following we use the shorthand dw = dx dt dy ds. Note that 
V '(u(x, t), v(y, a)) £>'>(-,*) - v(; S )](x, y) < &' r [ V (u(; t),v(; a))](x, y). 

Moreover, using the change of variables (x, y) — » (x — z,y — z), 

[[ [[ i>{x,y,t,a)&> r [r 1 {u{;t),v{;3))]{x,y)dw 

J J D T JJ D t 

= / / / / ri(u(x,t),v(y,s))il>(x- z,y- z,t,s) dwdfj,(z) 

J\z\>rJo Jz+AJO J z+A 

// // r](u(x,t),v(y,8))tp(x,y,t,a) dwdfi(z), 

J\z\>r JJ D T JJ D t 

which by periodicity and the definition of C*-^' r equals to 

/ / / / / V(u(x,t),v(y,a))ip(x - z,y- z,t,a) dwdfi(z) 

J\z\>rJo JaJo JA 

// // V(u(x,t),v(y,a)) ip(x,y,t,s) dwdfi(z) 

\z\>r JJ D T JJ D t 

r)(u(x, t), v(y, a)) £*'"' r [V(-, ■, t, a)](x, y) dw. 



1 D T JJ D T 

Therefore we have proved so far that 



r}(u(x, t), v(y, a)) (d t + d s )ip(x, y, t, s) 

' D T JJ D T 

+ q{u{x, t), v(y, a)) ■ (d x + d y )ip(x, y, t, s) 
+ V{u{x, t), v(y, a)) C*^ [ip{-,y, t, a)] (x) 
+ V {u{x, t), v{y, a)) Cp^(x, •, t, a)](y) 
+ V{u{x, t), v{y, a)) £*^ r [ip(-, •, t, a)](x, y) 
+ V(u(x, t), v(y, a)) 7^ • (d x + d y )ip(x, y, t, s) dw > 0. 
We now send r — > 0, remembering the definition of 7^ and defining 

C*^(;-)](x,y) 



J 

J\z 



}(x -z,y- z)- 4>(x, y) + z ■ (d x + d y )<j>{x, y) 1\ z \<i d/j,(z). 

l\z\>0 

The result is 



n{u{x, t), v(y, a)) {d t + d a )ip(x, y, t, s) 

' D T JJ D T 

" ■ ' ■ • 1 ' + q(u(x, t),v(y, a)) ■ (d x + d y )ip(x, y, t, s) 

+ 7]{u{x,t),v(y,a))jC*^(;;t,a)}(x,y) dw > 0. 
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To conclude, we show how to derive the /^-contraction (I2.2j) from this inequality 
by choosing the test function tp as 

(A.2) M Xi y it ,8) = & p U^^w s (j^)<l>(t), P,S>0, 

where lus(t) = for a nonnegative u> £ C^°(R) satisfying 

u>(-t) = w(r), w(r) = for all |r| > 1, and J uj(t) At = 1, 

while uj p (x) = ujp{x\) • • ■ uj p (x<i) with uj p {-) such that 

u) p (r) = w,,(t + 27rA;) . 

Note that w p is periodic in each coordinate direction. By a direct computation, 

{d t + d s )i>{x,y,t,s)=w p fcjMwa (j^-J 4>\t), 
(d x +d v )ip(x,y,t,s) = 0, 
£*[1>(;;t,a)](x,y)=0. 
Thus, with this test function tp at hand, inequality (|A.1|) becomes 

(A.3) \u(x, t) - v(y, s)\u, p (^^j co s cf>'{t) Aw > 0. 

We then go to the limit as (p, 8) — > to find that 

(A.4) // \u(x,t)-v(x,t)\(j)'(t) AxAt>0. 

J J D T 

To conclude the proof we now take <fi = Xm f° r 

(A.5) *„(*)=/ K( r - *i) - "/.(t - *a)) dr, < t x < t 2 < T. 

Loosely speaking, the function ^ is a smooth approximation of the indicator func- 
tion l(t!,t 2 ) which is zero near t = and t — T for /i > small. Since 

^(t) = w„(t-ti) -£«;„(* -t 2 ), 
inequality (IA.4I) reduces to 



\u(x, t) — v{x, t)\ ujfiit — t 2 ) dt dx < // |u(a;,t) — i>(a;,i)| ui^it — t\) At Ax. 

JJ Q T 

By the integrability of u and v and Fubini's theorem, the function 

= / |n(ar,t) dx € ^(0,7), 

and we may write the above inequality as a convolution 

$*w M (i 2 ) < <£>*oj p (t 1 ). 

By standard properties of convolutions, $ * — > a.e. i as fi — > 0. Hence, 

\\{u-v)(-,t 2 )\\ L i {A} < ||(it-u)(-,ti)||z,i(A) for a.e. t u t 2 £ (0,T). 

Finally, the theorem follows from renaming t 2 and using part Hi) in Definition 12.11 
to send t\ — > 0. 
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Appendix B. Proof of Theorem 12.41 

The vanishing viscosity problem (|2.3[) has a unique classical solution u e for s > 0, 
see Remark l2.6l If we multiply (|2.3[) by rj'(u e ) for any smooth convex function 77, use 
standard manipulations on the conservation law part combined with the inequalities 



r/(u e )£"[u t ] = r/K) ^[u e ] + C»> r [u e ] \ < C^{u e )] + r)'(u e )£»> r [u e ], 

r]'(u e ) Au e = Ar]{u £ ) -er)"(u e )\d x u e \ 2 < Arj(u e ), 

we find, after integration against any nonnegative test function 0, that u e satisfies 
the (entropy) inequality 

■q(u e7 k)d t ip + q(u e ,k) ■ d x ip + rj(u e , k) C*^[ip] +rf(u e ,k) £ M ' r [u £ ] 

D T 

+ r](u e , k) 7^ • d x (f + e rj{u e , k) Aip dx dt > 0. 

From this inequality we proceed as in the proof of the Z/^-contraction (Theorem 
I2.2j) . We take u = u(x, t), u e = u e (y, s), and find the inequalities 



rj(u(x, t),u e (y, s)) d t ijj(x, y, t, s) 

! D T JJ D t 

+ q(u(x, t), u e (y, s)) ■ d x ip(x, y, t, s) 

+ v( u ( x > t),u e {y, s )) ^bPi-, y, t, s )K x ) 

+ T)'(u(x, t), u e (y, s)) C^' r [u{-, t)](x) ip(x, y, t, s) 
+ r)(u(x, t), u e (y, s)) 7^ • d x ip(x, y, t, s) dw > 0. 

and 

// // ri{u{x,t) 1 u <L {y,s))d s ^{x,y,t,s) 

JJ D T JJ D T 

+ q(u(x, t), u e (y, s)) ■ d y ij}(x, y, t, s) 

+ v(u(x, t), u e (y, s)) £*^[ip(x, ■, t, s)](y) 

- v'(u(x, t), u e (y, s)) C^ r [u e (-, s)](y) ip(x, y, t, s) 

+ r)(u(x,t),u e (y,s))jj l ■ d y ip(x,y,t,s) dw 

+ erj(u(x,t), u e (y, s)) A y ^(x, y,t, s) dw > 0. 

As in the proof of Theorem 12.21 we add and manipulate these to get (see (jA.l 

rj(u(x, t),u e (y, s)) (d t + d s )ip(x, y, t, s) 



'D T JJ D t 

+ q(u(x, t),u e (y, s)) ■ (d x + d y )ip(x, y, t, s) 
+ r}{u{x, t),u e (y, a)) £*>1#, •, t, s)](x, y) 
+ er)(u(x,t),u e (y, s)) A y ip(x,y,t, s) dw > 0. 

We now take the test function i\) as in (TA~2|) and find that (see (|AT3|) ) 

(x — y \ ft, — ■ s \ 



(B.l) -JDtJJDt 



<e // ri(u(x,t),u e (y,s)) A y ip(x,y,t, s) dw. 

JJ D T JJ D t 
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After an integration by parts, the right-hand side (R.H.S.) is bounded by 
R.H.S. < e // // d v \u(x,t) - u e (y,s)\ d v i>{x,y,t,s) &w 

JJ D T JJ D T 

- e // // \ d v u ^V' s )\ IdyV'foj/jM)! Aw 

JJ D T JJ D t 

< CT \u \bv(A)-, 
P 

where the last inequality is a consequence of the estimate \u £ {-, 01 w(A) < |uo|w(A) 
and (|A.2p . 

To estimate the left hand side (L.H.S.) of (jB.ll) . note that 

- K(y,s) - u(x,t)\4>'(t) 

> -\u e (x,t)-u(x,t)\<j>'(t)-\u e {x,s)-u e (x,t)\\<p'(t)\ - \u e (y,s)-u e (x,s)\\<P'(t)\, 
and that 

\u e (x,s) - u e (x,t)\ui p (~Y^j u& (~~2~~) dW ^ ° 



I D T JJDr 

and 



(x — y \ ft — ^ \ 
~9~) U5 \2~ ) ^'(^ dw - CT \ u o\bv P- 

Hence we conclude after sending 8 — » that 

u e (x, t) - u(x, t)\ 4>'{t)dxdt ~Cp< L.H.S. (< R.H.S.). 



I D T 

The results then follows by setting p = yfk and (j) = \^ as m (|A.5[) . and conclude 
as in the proof of Theorem 12.21 Sending fi — > 0, setting t 2 — t, and using part Hi) 
in Definition 12. II to send t\ — > 0. 
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